#Read in Oil Spill Data
oil_spill <- read_sf(here("Oil_Spill_Incident_Tracking_%5Bds394%5D-shp"), layer = "Oil_Spill_Incident_Tracking_%5Bds394%5D") %>%
clean_names() %>%
rename(county = localecoun)
#Read in California County Boundaries
ca_county <- read_sf(here("ca_counties"), layer = "CA_Counties_TIGER2016"
) %>%
clean_names()
#Set county crs to same as oil spill
ca_county <- st_transform(ca_county, st_crs(oil_spill))
tmap_mode("view")
tm_shape(oil_spill) +
tm_dots("county")
tm_shape(ca_county) +
tm_fill("name", palette = "YlOrBr", show.legend = FALSE) +
tm_shape(oil_spill) +
tm_dots(aes(color = "brown"))